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Abstract: 

We report the variation of tiie spectral energy distribution (SED) of blazars as a function of source activity, based on 
available, simultaneous multi-wavelength (MWL) observations of BL Lac objects. We use a fully automatized min- 
imization procedure, instead of the commonly used eye-ball fit, to model the data sets with a one-zone Synchrotron- 
Self-Compton (SSC) model. The obtained SSC parameters are then analyzed as a function of source luminosity, and the 
correlation between parameters is shown. Possibilities of improving the present observational and modeling status of BL 
Lac objects are also discussed. 
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1 Introduction 



scattered by the same electron population responsible for 
the synchrotron emission ||2l . 



Active Galactic Nuclei (AGNs) are believed to be the most 
powerful energy sources in the universe. It is usually as- 
sumed that at the AGN center there is a super massive black 
hole (SMBH) surrounded by the accretion disk: a beamed 
emission perpendicular to the disk is expelled in opposite 
directions in ^10% of the cases (radio-loud AGNs). The 
energy required for the beamed emission might be com- 
ing from the rotating black hole and the fully ionized ac- 
cretion disk [1 1. AGNs are phenomenologically classified 
according to their viewing angle to the observer If the jet 
is pointing very close to the line of sight to the observer 
the source is called a blazar. Blazars include Flat Spectrum 
Radio Quasars (FSRQ) and BL Lac objects. 

Blazar emission is dominated by a non-thermal continuum. 
The spectral energy distribution (SED) shows two peaks. 
For the subclass of high-peaked BL-Lac objects (HBLs) 
the lower-energy peak is in the UV/X-ray range and the 
higher-energy peak is in multi-GeV bands; the so called 
low-peaked BL Lac objects (LBLs) have, instead, the low- 
energy peak in the IR/optical band and the high-energy 
peak in the multi-MeV band. The low-energy peak is due 
to synchrotron radiation. According to the leptonic model, 
the high-energy peak is produced by inverse Compton (IC) 
upscattering of lower-energy photons off relativistic elec- 
trons: in particular the so called synchrotron self Compton 
(SSC) model suggests that the synchrotron photons are up- 



^ minimized one-zone SSC model 
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In our analysis we will employ a specific one-zone SSC 
model |3|. This model assumes that the source is a spher- 
ical blob of plasma of radius R, moving with a Doppler 
factor 6 towards the observer at an angle 9 with respect to 
the line of sight, threaded with a uniform tangled magnetic 
field of strength B. The injected relativistic particle pop- 
ulation is described as a broken power-law spectrum with 
normalization K, extending from ^^in to 7max with indices 
Til and 712 below and above the break Lorentz factor 7brcak- 
So, the total number of free parameter of the model is 9: 6 
parameters specify the properties, and 3 the environment, 
of the emitting electrons. The one-zone SSC model can 
be fully constrained by using simultaneous multi-frequency 
broad-band data. 

In the present work we have kept 7inin fixed and equal to 1, 
as usually done in the literature. The determination of the 
remaining parameters has been performed by finding their 
best values and uncertainties from a -minimization in 
which multi-frequency experimental points have been fit- 
ted to the SSC model described in Minimization has 
been performed using the Levenberg-Marquardt method 
|l4l, which is an efficient standard for non-linear least- 
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Table 1: Log of SED data sets. 
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squares minimization that smoothly interpolates between 
two different minimization approaches, namely the inverse 
Hessian method and the steepest descent method. 

A crucial point in our implementation is that from |3 | we 
can only obtain a numerical approximation to the SSC 
spectrum, in the form of a sampled SED. On the other 
hand, at each step of the loop the calculation of re- 
quires the evaluation of the SED for all the observed fre- 
quencies. Although an observed point will likely not be 
one of the sampled points coming from |3 1, it will fall be- 
tween two sampled points, so that interpolation can be used 
to approximate the value of the SED. At the same time, 
the Levenberg-Marquardt method requires the calculation 
of the partial derivatives of with respect to the SSC pa- 
rameters. These derivatives have also been obtained nu- 
merically by evaluating the incremental ratio of the with 
respect to a sufficiently small, dynamically adjusted incre- 
ment of each parameter. 

3 Data selection 

In order to study the behavior of SSC parameters for each 
source in different emission states, we selected sources 
with at least 2 available simultaneous SED data and with 
known redshift. The MWL campaign period of each data 
sets used for this study is given in Table 1 . 

4 Results 

Our results suggest that in Mkn421, B decreases with 
source activity whereas 7break and 5 increase (Fig. 3 top). 
This can be interpreted in a framework in which the syn- 
chrotron power and peak frequency remain constant with 
varying source activity by decreasing B and increasing the 
number of low-energy electrons. This mechanism results 
in an increased electron-photon scattering efficiency and 
hence in an increased Compton power Other emission 
parameters appear uncorrected with source activity. In 




45 4a 

Log L [erg g-i] 




45 4B.B 4S 

Log L [erg s-i] 



Figure 3: The SSC fit parameters jbreak, S and B versus 
source luminosity in Mkn 42 1 

Fig. 3 (bottom), the -B-7break anti-correlation results from 
a roughly constant synchrotron peak frequency. The B-5 
correlation suggests that the Compton peak of Mkn 421 is 
always in the Thomson limit. The (5-7bicak correlation is an 
effect of the constant synchrotron and Compton frequen- 
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Figure 4: The correlation between the SSC fit parameters 
in Mkn421 



To test the goodness of the fits we have apphed the 
Kolmogorov-Smimov (KS) test for normality of the resid- 
uals to all the SED fits presented in this paper. A standard 
application of the test shows that the residuals are not nor- 
mally distributed in all cases, i.e. for most samples the KS 
test fails at 5% significance. This could happen because 
of two reasons, (i) The one zone SSC model, which con- 
tains two distinct physical processes (synchrotron emission 
and its Compton upscattered counterpart) may only be an 
approximation to the real process: additional subtle effects 
may also enter the modeling of blazar emission, (ii) In our 
blazar data sets the uncertainties associated with 7 data are 
much larger than those associated with the optical and X- 
ray data. Both these observations suggest that the test could 
give different results if applied separately at low- and high- 
energy data. We found that even though the KS test fails 
on the whole data sets, it instead confirms normaUty of the 
residuals at 5% confidence level when applied piecewise to 
low and high energy data. 

With better and more complete VHE data, it would be pos- 
sible to reduce the uncertainties in the parameters, to ob- 
tain a higher overall significance in the fit and, therefore, 
a clearer picture of the changes and properties and envi- 
ronment of the emitting electrons as a function of source 
activity. 
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cies of the radiation emitted by a plasma in bulk relativistic 
motion towards the observer. 

The analysis of the remaining sources is in progress. How- 
ever, it should be noted that simultaneous MWL campaigns 
for other blazars are fewer than for Mkn 421 . Hence we en- 
courage more MWL campaigns on blazars at their different 
emission level. 
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Figure 1: Mkn421 in different emission states. 
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Figure 2: Other BL Lac objects (listed in Table 1) in different emission states. 



